Soret Motion of a Charged Spherical Colloid 
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The thermophoretic motion of a charged spherical colloidal particle and its accompanying cloud of 
counterions and co-ions in a temperature gradient is studied theoretically. Using the Debye-Hiickel 
approximation, the Soret drift velocity of a weakly charged colloid is calculated analytically. For 
highly charged colloids, the nonlinear system of electrokinetic equations is solved numerically, and 
the effects of high surface potential, dielectrophoresis, and convection are examined. Our results are 
in good agreement with some of the recent experiments on highly charged colloids without using 
adjustable parameters. 
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Introduction. A temperature gradient applied to a 
fluid mixture causes relative transformation of its com- 
ponents; some condensing in the hotter and some in 
the colder side. This phenomenon, known as the Soret 
effect or thermophoresis, has been studied for nearly 
150 years [1 1 and has been observed in a variety of 
systems [a~IE H H @, 0. While the existence of 
such a response can be well formulated in the frame- 
work of non-equilibrium thermodynamics Q, in many 
cases its microscopic nature remains poorly understood 
E3, El, E El 13 ■ The manifestation of this effect in 
charged colloidal solutions is particularly puzzling: while 
it is usually observed that colloids condense in the colder 
side, experiments with opposite results also exist Q and 
it appears that the tendency could change with the vari- 
ation of salt concentration and temperature [H, Q| ■ 

In a pioneering work in 1981, Ruckenstein suggested a 
model for Soret motion of a single charged colloidal par- 
ticle [9j. Using Boltzmann distribution for ion densities 
and the Debye-Hiickel approximation for the electrostatic 
potential, he was able to solve the hydrodynamics equa- 
tions and find the drift velocity of a colloid in a tempera- 
ture gradient. The 2002 experiment by Piazza and Guar- 
ino on Sodium Dodecyl Sulfate (SDS) micelles solutions 
provided an opportunity for verification of this model Q • 
The SDS micelles in the experiment were highly charged, 
as they had an estimated saturated structural charge of 
Z ~ 75 (in electronic unit) [3 and a radius of a = 2.5 
nm, and the thickness of their double-layer — as set by 
the Debye length — varied between 0.4 nm to 2.5 nm [3j. 
The experiment was therefore outside the region of va- 
lidity of the Ruckenstein's theory, which was restricted 
to weakly charged colloids with thin double-layers. Pi- 
azza and Guarino showed that it was possible to get a 
reasonable fit to the experimental data for the depen- 
dence of the Soret coefficient on the salt concentration 
using Ruckenstein's formula, provided they assumed an 
increased radius of a = 3.5 nm and a reduced charge 
of Z = 17 [3]. The apparent smaller charge of the col- 



loid could be interpreted as the renormalized charge as 
described by Alexander et al. [3- However, this attri- 
bution is not entirely justified as the concept of effective 
charge is defined through the asymptotic form of the elec- 
trostatic potential of the colloid far away from its surface 
[3 - while Ruckenstein's model deals with the electric 
field inside the thin double-layer. In fact, the value of 
the effective charge Z = 17 used by Piazza and Guar- 
ino are obtained from measurements on inter-colloidal in- 
teractions [3, and electrophoretic observations on SDS 
micelles does not seem to confirm this value 17J. More 



recently, Putnam et al. have performed experiments on 
highly charged T4 lysozymes and found results that do 
not seem to be explained satisfactorily by any of the ex- 
isting theories on single-colloid thermophoresis 0]. An- 
other recent experiment by Duhr and Braun [B[ probed 
charged particles with low surface potentials, and found 
a new scaling relation for Soret motion of charged col- 
loid which docs not agree with Ruckenstein's prediction. 
They showed that it is possible to explain their obser- 
vations, using the Gibbs enthalpy of the charged colloid 
Q (or the irreversible work needed to construct it [3)- 
The controversy was later clarified to some extent but the 
work of Astumian [3 > who suggested that while Rucken- 
stein's model deals with the deterministic steady motion 
of a charged colloid, the theory by Duhr and Braun ad- 
dresses its fluctuation-induced stochastic motion, which 
is a separate contribution, and a complete picture should 
involve both of these aspects simultaneously. 

Here, we focus on the deterministic motion of a charged 
spherical particle in a temperature gradient, and consider 
both weakly charged and highly charged cases taking into 
account effects such as the temperature dependence of 
solvent electric permittivity, convection, and nonequilib- 
rium co/counterion redistribution. We examine the be- 
havior of system using Debye-Hiickel approximation for 
weakly charged colloids, and provide an analytical result 
for colloid drift velocity with arbitrary double-layer thick- 
nesses. For highly charged colloids, we solve the nonlin- 
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ear set of coupled equations numerically, and find that 
the Soret coefficient has a non-monotonic dependence on 
the surface (zeta) potential of the colloids, an effect that 
has been indeed observed by Putnam et al. Q ■ We also 
examine the relative importance of various contributions 
in the different regimes. For the experiment with SDS 
micelles 0], we show that our systematic approach can 
yield results in reasonable agreement with the experi- 
mental data, using a realistic bare value for the micellar 
charge. 

The model. We consider an ionic solution of different 
species with concentration Ci and valence qi, which cre- 
ate an electrostatic potential cf> through Poisson's equa- 
tion —V ■ eV<fi = 4-7T YJi QiCi. This solution is subject to 
a gradient in temperature which causes e to change in 
space. We consider weakly varying temperature fields, 
so use linear response theory to study the deformation 
of the double-layer and solvent flow. The Soret motion 
of a charged colloid is a steady motion which means that 
in the colloid framework, we should look for a stationary 
solution for the fluid velocity field V. This fluid veloc- 
ity is governed by the Navier-Stokes equation and the 
incomprcssibility constraint V • V — 0. In the limit of 
low Reynolds number, this yields the stationary Stokes 
equation: — -q\7 2 V — —VP — V</>Es + / Here rj is 
the viscosity of the solvent, and P is its pressure. The 
body force / is the dielectrophoretic force, which comes 
from the net force experienced by electric dipoles of wa- 
ter molecules because of spatial variation in the electric 

field [H3]. It reads /= V (^-|£| 2 ) - ^\E\ 2 Ve 

where E = — V0. The first term on this expression is a 
complete derivative and can be absorbed in the fluid pres- 
sure (i.e. P — > P — ^!-|i?| 2 ). Assuming that the solvent 
is incompressible (i.e. p = po), changes in its permit- 
tivity e(po,T) is due to the temperature gradient only: 
—±\E\ 2 Ve = a^\E\ 2 ^- where a = -d\ne/d\nT is 
about 1.35 for water in room temperature [2p| . Finally, 
each ion species is subject to a conservation law V- J, = 0, 
with ionic current density: 

Ji = -DiVQ - t MC iqi Vcj) + VQ - DidS^VT (1) 

where pi is the ith type of ions mobility and Di = piksT 
is its diffusivity. 5™ is the co/counter-ions Soret coef- 
ficient. Here we only focus on 1 : I electrolytes, so we 
assume that coions and counterions have equal Soret co- 
efficients [in. 

In the absence of any temperature gradient and 
fluid motion, Eq. JI} is simplified to its first 
two terms. Then, Ci = Co exp[— q^/ksT] yields 
Ji = 0, and satisfies V • Ji = 0. In the pres- 
ence of a temperature gradient, however, not only 
we have to consider all of Eq. (p} terms, but also 
note that a Boltzmann weight form for G no longer 
makes — DjVCj — piCiqiS '4> vanish; instead, it yields 
—DiCi(qi4>/kBT 2 )VT. We suggest to use the following 
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measures the deviation of the concentration 
from a Boltzmann weight form, and it contains contribu- 
tion from the convective term in Eq. {T]) as well as the 
aforementioned term of —DiCi(qi4>/kBT 2 )VT. 

Our aim is to find the colloid drift velocity, so we focus 
on the Stokes equation that governs fluid velocity. The 
electric and dielectrophoretic forces on the right hand 
side are acting as source terms that induce fluid motion. 
To first order in temperature changes, the source terms 
simplify to 
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where Cf = Co exp[— qi4>o/kBTo] and Eq are the non- 
perturbed values of ion density and electric field. Here, 
we have also extracted another complete derivative 
(i.e.— V V\ qjCf64>), which can be absorbed in the fluid 
pressure term jl 81 ] . 

Weakly Charged Colloids. For q<p/k^T <C 1, we can 
use the Debye-Hiickel approximation and solve this sys- 
tem of equations analytically The drift velocity of 
a colloid with radius a is found as 



Vdrift = 



(3) 



48tt t]T 

x { (I + T S£ n )F( K a) - G(ko) + a [2 - F(ko)} } , 

where <ps = Zq/ea(\ + na) is the zeta potential of 
the surface of the colloid and k = y / 8nq 2 Co/ekBTo 
is the inverse Debye length [22||. In this equation, 
F(x) = 2x- 4x 2 e 2x E 1 (2x) and G(x) = f [x(l + x)(12 - 
x 2 )e x Ei(x) + 8~llx + x 3 ~24xe 2x E 1 (2x)}, with Ei(x) = 
f°° e~ s ds/s. Equation ([3]) is arranged in the form of 
Ruckenstein's result for colloid drift velocity 0] multi- 
plied by a correcting factor. Each term in this correcting 
factor corresponds to one of the source terms in Eq. (|2|) 
(i.e. a [2 — F(na)] corresponds to the dielectrophoretic 
term, -G(na) to n t ST/T term, and T S^ n F( K a) to 
ions Soret term). The remaining F(ko) corresponds to 

the -V0o Ei <li C i THr a ST / T o) term, which corrects 
Ruckenstein's result [j| for arbitrary double-layer thick- 
nesses. 

The limit of a thick double-layer (i.e. na -C 1) or low 
ionic strength is particularly interesting, as for na = 
we have F(na) = G(na) = [HI], and the only contri- 
bution in the correcting factor of Eq. (J3j) will be 2a. It 
means that with low ionic density, the dielectrophoretic 
force that is the force on water molecules [HI, plays the 
dominant role in the phenomenon. In addition, Rucken- 
stein's formula will be still applicable, if we multiply it 
with a constant 2a, which can be presented in terms of 
a renormalized surface potential (f) s — \[2a<^s or charge 
Z = V2a~Z [H. 
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FIG. 1: (color online), (a) TSt versus co for a colloid with 
radius a = 2.5nm, but different charges. Mean temperature 
is T — 25° C, and we extrapolated existing data [2~l| to obtain 
(Na+ and CL") ions Soret coefficient: = 0.99 - 2 x f0~ 3 
K _1 (depending on salt density), (b) Ratios of the dielec- 
trophoretic term (a-term) to the .F-term (dashed blue line) 
and G-term/F-term (solid red line), versus \q<j>s /kBT\, for 
kb = 2. \q<f>s/k B T\ = 3.74 corresponds to the Z = 50 - 60, 
our suggested charge for a micelle in the experiment of Ref. 
Q|. The dashed-dotted (orange) curve corresponds to when 
the convection term is artificially switched off. 



Highly Charged Colloids. The Dcbyc-Hiickcl approx- 
imation is not valid for highly charged colloids and one 
needs to fully take account of the nonlinearity of the 
electrostatics. To this end, we have solved the above 
governing equations numerically. Figure [TJi shows the 
Soret coefficient TSt ~ Klrift a s a function of na for 
various values of structural charge, corresponding to 
the SDS micelles used in the experiment of Ref. Q. 
We find that while the Soret coefficient initially in- 
creases with the structural charge — in agreement with 
the Dcbyc-Hiickel limit behavior where it is proportional 
to <j> 2 s ~ Z 2 /{1 + na) 2 — this trend is reversed at suffi- 
ciently high values of the charge and the Soret coefficient 
starts to decrease until it changes sign and becomes neg- 
ative. We can understand this behavior better, if we 
track the contributions of the different source terms in 
Eq. ([2]) to the Soret coefficient. Following our analytical 
approach [that led to Eq. (|3j)]. we name the contribution 
of -V0o E< g i C l o (g i 0o/fcsTo)(5T/T o ) term in Eq. © as 
the F-term, that of the dielectrophoretic term as the a- 
term, and finally that of fii's as the —G-term (note minus 
sign of G(na) in Eq. ([3])). Figure [TJd shows the relative 
importance of these contributions as a function the sur- 
face potential. While for weakly charged colloids (i.e. 
\q<ps /k-BT\ <C f), the a-term is the dominant contribu- 
tion, for \q(j)s /kBT\ > 1, the situation changes in favor 



of the G-term, which eventually leads to negative Soret 
coefficient at high surface potentials, interestingly, for 
the surface potentials attributed to the SDS micelles of 
Ref. 0, the two terms balance each other. 

A similar non-monotonic behavior has recently been 
observed by Putnam et al. for lysozyme proteins with 
various surface charges/potentials Q- Figure [2jt shows a 
comparison between their observed TSt versus lysozyme 
surface potential, and the prediction of our theory with- 
out any adjustable parameters. In this calculation, we 
have ignored the contribution of the Soret coefficient 
of the ions, as we were not able to obtain relevant ex- 
perimental values for these quantities. The encourag- 
ing agreement suggests that for highly charged colloids 
it will not be sufficient to treat the double layer close 
to the surface within the equilibrium Poisson-Boltzmann 
formulation, and it is necessary to fully take account of 
the nonlinear nonequilibrium effects in the profile of the 
ion cloud around the colloid. 

It is also interesting to probe the role of convection 
in our results. The convection term in Eq. (QJ affects 
fii's and consequently the G-term. The dashed-dotted 
(orange) curve in Fig. [TJd shows the ratio G-term/F- 
term, with convection artificially turned off. We ob- 
serve that while for \q<ps /ksT\ <C f the convective term 
has a small effect (proportional to |(?0s/fcsT| 2 [HI]), for 
higher surface potentials the convective motion will be 
non-negligible, but always finite. 

We also do numerical calculation for experiment of Pi- 
azza and Guarino, and use electrophoresis data [f7| to 
find micelles charge [H| . The result is presented in Fig. 
[2b, which only shows an agreement with the experiment 
for Kfl > 3. Since we have established (see Fig. [T]) that 
changing the structural charge does not lead to a bet- 
ter harmony, we can conclude that effects other than the 
ones considered here systematically should control the 
behavior of the system for na < 3. 

Discussion. The experiment performed by Duhr and 
Braun unlike most other experiments, have been done 
for colloids within Dcbyc-Hiickel regime. Therefore, it is 
natural to expect that our analytical results should agree 
with its findings. Figure [2fc compares our result (based 
on Eq. Q) with its data, and does not show a satis- 
factory agreement. We have examined possible correc- 
tions coming from fluid slip velocity on the colloid sur- 
face and mismatch in the thermal conductivities of the 
colloid and the fluid, and found that they cannot im- 
prove the situation 12411 . To resolve this discrepancy it 
has been suggested [13J that one should differentiate be- 
tween the regime where hydrodynamic or deterministic 
components of the motion are dominant, and the regime 
in which the system feels a local thermal equilibrium and 
it is the stochastic motion which is the dominant player. 
It is expected that the behavior of the Soret coefficient as 
a function of various parameters such as the radius, tem- 
perature, etc. will be significantly different in these two 
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FIG. 2: (color online). (a) Comparison of Putnam et 
al. [3] data (red) for mutant T4L lysozymes with Z = 
+3, +5, +7, +9, and Rh ~ 1.8nm, with our numerical pre- 
diction (dashed-dotted green curve), (b) Comparison of the 
experimental data of Piazza and Guarino [j| with our theoret- 
ical prediction. Depending on salt concentration, we assume 
a varying micelle charge of Z = 50 — 60 [l8|[. a varying ions 
Soret Coefficient S£ n = 0.99-2 x 10" 3 K _1 ^nd a fixed ionic 
diffusion coefficient D = 1.33 x 10 9 nm 2 /s [24|]. Our data is 
lifted by +6.547, to obtain the best fit with four last points, 
(c) Comparison of our analytic theory with Duhr and Braun 
data [H for polystyrene spheres with different radii. 

regimes. This means that a theoretical formulation that 
can account for experiments corresponding to the hydro- 
dynamic or deterministic regime is almost bound to fail 
to account for experiments in the stochastic regime as 
they would be mutually exclusive, a view which is also 
supported by recent experimental evidence [25|. Further 
work is needed to fully clarify this picture, ideally in the 
form of a theoretical formulation that includes both ef- 
fects and can show the crossover from one behavior to 
the other in a systematic way. 

In conclusion, we have presented a systematic theoret- 
ical analysis of the Soret motion of charged colloids both 
in the weakly and strongly charged regimes. Using con- 
tributions from different physical sources, we have found 
competing tendencies that are entirely due to the non- 
linearity of electrostatics and the noncquilibrium redis- 
tribution of the ion clouds around the colloid. 

We are grateful to E. Bringuier, R. Ejtchadi, F. Miri, 
and P. Pincus for helpful discussions; and A. Naji, S. A. 
Putnam, and G. Wong for their assistance. 
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